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EFFECTS OF TEMPERATURE DISTRIBUTION AND ELASTIC PROPERTIES 
OF xMATERIALS ON GAS-TURBINE-DISK STRESSES 

By Ahthcs G. Holms and Richard D. Faldetta 


SUMMARY 

Calculations were made to determine the influence of changes 
in temperature distribution and in elastic material properties on 
calculated elastic stresses for a typical gas-turbine disk. 

Severe temperature gradients caused thermal stresses of 
siiflicient magnitude to reduce the operating safety of the disk. 
Small temperature gradients were found to be desirable because 
they produced thermal stresses that subtracted from the centrifu- 
gal stresses in the region of the rim. The thermal gradients 
produced a tendency for a severe stress condition to exist near 
the rim but this stress condition could be shifted away from the 
region of blade attachment by altering the temperature distri- 
bution. 

The investigation of elastic material properties showed that 
Centrifugal stresses are slightly affected by changes in modulus 
of elasticity, but that thermal stresses are approximately pro- 
portional to modulus of elasticity and to coefficient of thermal 
expansion. Where thermal stresses are significant, use of 
materials with low moduli of elasticity and low coefficients of 
expansion would be desirable. Design calculations for disks 
to be operated at elevated temperatures would require accurate 
data on the values of these properties at all temperatures en- 
countered in order that accurate calculations could be made. 
Centrifugal tangential stresses at the rim were decreased by 
increases in Poisson’s ratio but the effect of changes in Poisson’s 
ratio on the total stresses was small. Elevated-temperature 
measurements of Poisson’s ratio would not, therefore, be 
e.ssential for disk-stress calculations. 

INTRODUCTION 

Gas-turbine disks are normally operated at such high 
temperatures that the materials used are at a low-strength 
level. The hot gases contact the blades and the rim of the 
turbine rotor and thus maintain the rim at a high tempera- 
ture. Various cooling methods have been used to reduce 
the temperature of the disk, but because the rim is always 
in contact with hot gases, it remains at a high temperature, 
whereas cooling decreases the temperature of the central por- 
tion of the rotor and thus increases the temperature gradients. 
These gradients are sources of thermal stresses that cause 
the stress distributions in gas-turbine disks to differ wideK 
from those encountered in steam turbines, as steam turbines 
are mainly subjected to centrifugal stresses with small 
temperature gradients. An analysis of the types of stress 


to be e.xpected in aircraft gas turbines requires a considera- 
tion both of the centrifugal stresses and of the stresses 
resulting from temperature gradients. 

In order to investigate the net effects of superposed thermal 
and centrifugal stresses, a number of arbitrarily selected 
cooling conditions were assumed to exist in a typical disk 
and the corresponding thermal and centrifugal stresses were 
calculated. The margin of safety was deter min ed by 
algebraically adding the centrifugal and thermal stresses 
present at any point and comparing the resulting stress state 
with the rupture strength of the material at the temperature 
of the point. This comparison of the strength of the material 
with stress conditions was considered to afford an indication 
of the relative desirability of various temperature dis- 
tributions. 

Elastic properties of materials such as Poisson’s ratio, 
modulus of elasticity, and coefEcient of thermal expansion 
vary with temperature. Aircraft gas turbines are usually 
operated with a radial variation of temperature in the disk. 
Accurate stress calculations would therefore require accurate 
data on the elastic properties of materials at all operating 
temperatures and a method of stress calculations woidd have 
to be used that would account for the point-to-point varia- 
tion of material properties in a disk. 

In order to determine which properties significantly affect 
the stresses in gas-turbine disks and to determine the t 3 ^pes 
of effect that result from varying these properties, calcula- 
tions of centrifugal and thermal stresses were made for a 
wide variation of values of elastic constants in the fall of 
1946 at the NACA Cleveland laboratorj^. These calcula- 
tions were facilitated by a method (reference 1) that is 
especially adapted to disks in which there are point-to-point 
variations of profile, temperature, and material properties. 

.ASSUMPTIONS AND METHOD OF ANALY'SIS 

Centrifugal stresses were calculated for a rotative speed of 
11,000 rpm. The disk profile used in all calculations is 
shown in figures that are the results of stress computations. 
In the calculations of the influence of changes in radial 
temperature distribution, the assiunption was made that no 
axial stresses or temperature gradients are present. (Axial 
stresses and temperature gradients would be expected in 
disks that are cooled on one side.) 

It was assumed that the radial temperature distribution 
could be represented a function of the form 
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where 

T temperature at radius r 

a, b, % constants determined by assumed temperature dis- 
tribution 
r radius 

This functional relation between temperature and radius was 
chosen as a convenient method of representing anticipated 
temperature distributions. Temperature measurements 
could prove other types of function to be more appropriate 
to particular turbines. 

Centrifugal stresses were calculated on the assumption 
that the mean radial stress at the rim due to blade loading 
is a tensile stress of 8500 pounds per square inch. This value 
was obtained by dividing the total centrifugal force at the 
root of the blades by the total rim periplnral area. Tem- 
perature distribution in the blades was assumed to have no 
influence on the disk thermal stresses and the radial thermal 
stress was assumed to be zero at the rim. 

The method of stress calculation used is given in reference 1 . 
It is a finite-difference method that permits calculation 
of the stresses present in disks in which temperature, thick- 
ness, modulus of elasticity, Poisson’s ratio, and coefficient of 
expansion are varied in an arbitraiy manner. In order to 
calculate only centrifugal stresses, the coefficient of thei'iual 
expansion was set equal to zero; in order to calculate only 
thermal stresses, the angular velocity was set equal to zero. 
The total stresses were found by algebraically adding the 
centrifugal and thermal stresses. 

Because aircraft power plants are required to operate 
under severe conditions for only short periods of time (such 
as during starting and take-off'), the 2-hour rupture strength 
was arbitrarily chosen as a criterion of the capacity of a 
material to withstand the stresses of aircraft gas-turbine 
disks. The 2-hour rupture strength is the value of uniaxial 
tensile stress that will cause failure of a tensile test specimen 
in 2 hours at a definite temperature. The high-temperatirre 
portion of figure 1 is a plot of 2-hour stress-rupture values 
taken from figure 46 of reference 2; the low-temperature 
portion of the figm'e was obtained by passing a smooth 
curve tlu'ough the elevated-temperatiu'e points and a room- 
temperature point (107,000 Ib/sq in.), which was taken from 
the tensile-strength curve of figure 44 of reference 2. The 
disk stresses are biaxial; hence, comparison of the severity 
of stresses at a point in a disk with the stress-rupture value for 





(a) No uniform tomporature region at rim. (b) Small uniform temperatui'O region at rim 
(o) Large uniform temporatui-eregion at rim. (d) Variable center tomporaluj'o. 

Figure 2.— Assumed radial temperature distribution relative to turbine-disk pronic os given 

by F^ar'-H*. 


the corresponding temperature requires the use of a single 
quantity that represents the biaxial stress state. For this 
purpose, an “equivalent tensile stress,” which is defined by 
the derivation given in the appendix, was calculated from 
the maximum-shear-strain-energy theory of faihu-e. The 
excess of a stress-rupture value over the equivalent tensile 
stress is called the margin of safety. This quantity, which 
was calculated on the assumption that the disk remains 
perfectly elastic, provides a means of comparing the stress 
conditions at various points of the disk and also provides a 
qualitative means of studying the net effect of changes m 
temperature distribution. Calculation of the margin of 
safety using elastic stresses is not a method of estimating 
failure conditions because a disk, which is locally overstressed, 
win yield plastically and so relieve the high stresses caused 
by stress concentration. 

EFFECTS OF TEMPERATURE DISTRIBUTION 

In order to investigate the manner in which stresses depend 
on temperature distribution, a variety of arbitrarily selected 
temperature distributions (fig. 2) was assumed and the 
corresponding stresses were calculated. 

TYPICAL STRESS DISTRIBUTION 


A qualitative presentation of the radial variation of com- 
ponents of centrifugal and thermal stresses, total stresses, 
equivalent tensile stresses, 2-hour rupture strength, and 
margin of safety is shown for a typical disk in figure 3. 
The data for this illustration of the manner in which the 


Figure l.— Two-hour stresj-rupture curve for solution-quenched 16-25-6 alloy, (Data from 

reference 2.) 
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(a) Centrifugal and thermal stresses. (b) Total stresses. 

(c) Two-hour stress-rupture curve and Cd) Margin o£ safety, 

equivalent tensile stresses. 


Fiotre 3,— Typical stre.ss distribution for rroBle shown. Speed, 11,000 rpm; bucket loading 
at rim, 85(X} pounds per square inch; center temperature, 600* F; rim temperature, 1200® F; 
tt^mperature distrCbution given by T*ar»-{-6; 


center; In the event that these stresses become excessive, a 
sudden rupture of the rotor would probably occur. Disl^ 
are therefore usually designed to have a larger margin of 
safety at the center than at the rim. 

EFFECT OF CHANGING TEMFEBATURE DISTRIBUTION 

The margins of safety that correspond to the several 
temperature distributions of figure 2 (a) are shown in figure 4. 
The margin of safety at the center is large and increases as 
n is increased but, as the temperature gradients in the 
region of the rim are increased by increasing n, the margin 
of safety at the rim is reduced. Figure 5 shows curves of 
the equivalent tensile stresses and of the margins of safety that 
correspond to the temperature distributions of figure 2 (d). 




Figitre I - Distribution of margin of safety for profile shown. Speed, 11,000 rpin; blade 
loading at rim, 8500 pounds per square inch; center temperature, C00° F; rim temperature, 
1200° F; temperature distribution giren by T=or"F*. 


components of stress combine were calculated for the temper- 
ature distribution of figure 2 (a) with n equal to 2. The 
convention of considering tensile stresses positive and com- 
pressive stresses negative was followed. Figure 3 (a) shows 
that the centrifugal stresses are all tensile but the thermal 
tangential stress takes on large compressive values. Because 
the tangential thermal and centrifugal stresses at the rim 
are of opposite sign, the total tangential stress at the rim is 
smaller than the thermal stress alone. The two stress sys- 
tems have been added in figure 3 (b) and it is seen that a net 
tangential compressive stress exists in the region of the rim. 
The equivalent tensile stress for the centrifugal stresses and 
for the total-stress system are shown in figure 3 (c). The 
2-hour stress-rupture values corresponding to the assumed 
temperatures are also shown. Figure 3 (d) shows the cor- 
responding margin of safety. The low margin, of safety at 
the rim is due to the presence of resultant compressive 
stresses. Excessive values of such stresses often cause 
plastic flow of the disk rim during engine operation and when 
the disk cools a system of tensUe stresses is set up that can 
cause rim cracking. Because such cracks usually progress 
rather slowly, serious trouble from this source can, in most 
cases, be forestalled by removal of the wheel from service. 
The thermal and centrifugal stresses are both tensile at the 




(a) Equivalent tensile stress. 

(b) Margin of safety. 

Figure 5. — Tariation in equivalent tensile stress and in margin of safety for profile shown hs 
temperature distribution is varied according to figure 2 (d). Speed, 11,000 rpm: blade loa<l- 
Ing at rim, 85(X) pounds per square Inch; rim temperature, 1200® F; temperature dlstribatfoa 
given by n*2. 
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In general, the reduced thermal stresses decrease the 
equivalent tensile stresses as the temperature of the disk 
center is raised ; however, the stress at the rim is higher for a 
center temperature of 1200° F than for a center temperature 
of 1000° F. The tangential centrifugal and thermal stresses 
are of opposite sign and, for the temperature distribution 
corresponding to a center temperature of 1200° F, there is no 
thermal stress to reduce the effect of the centrifugal stress; 
therefore, the equivalent tensile stress at the rim for a center 
tempera tiire of 1000° F is lower than that for a center tempei- 
ature of 1200° F. Figure 5 (b) shows a different order of 
desirability of center temperatures from that which would 
be exjiectcd from figure 5 (a) when the center of the wheel 
is considered. This difference exists because low center 
temperatures raise the material strength at the center and 
thus influence the margin of safety. Those results show that 
the margin of safety at any point is affected by temperature 
distribution in two ways; that is, the thermal stresses and 
the strength of the material are affected by changes in 
tempeniture. 

CONTROL OF STRESS DISTRIBUTION BY ADJUSTMENT OF TEMPER ATURE 

DISTRIBUTION 

The attachment of blades to the disk rim usually increases 
the local stresses. Centrifugal forces and thermal gradients 
that produce low stresses hi the region of the rim are therefore 
desirable. If the turbine is .so designed that hot gases are in 
contact witii a large part of the rim, a uniform temperature 




(a) Tangential thermal stre.sses. 

(b) Margin of safety. 

Figcre 6.— Distribution of tangential thermal stress and margin of safety for several tem- 
perature distributions of profile shown. Speed, 11,000 rpm; blade loading at rim, 8500 
pounds per square inch; center temperature, 600“^ F; rim temperature, 1200° F; temperature 
distribution given by 



(a) Various slopes. 

(b) Zero slope. 

(c) Constant, typical slope. 


Figure 7.— Assumed relations between modulus of elasticity and temperatm e for calculations 
of thermal and centrifugal stresses. 


region at the rim would be expected, as is illustrated by 
figures 2 (b) and 2 (c). In figure 6 (a), wdiieh shows the 
effect of a constant temperature near the rim, are plotted the 
tangential thermal stresses that accompany the temperature 
distributions of figures 2 (a), 2 (b), and 2 (c) with n equal to 
2. Figure 6 (b), which is a plot of the corresponding margins 
of safety when all thermal and centrifugal stresses are con- 
sidered, shows that the position of minimum margin of 
safety can he shifted by changing the temperature distribu- 
tion. One method of changing the temperature distribution 
would be to relocate the seal between the cooling air and the 
hot gases. Another method would be to alter the supiily of 
cooling air. In this manner, shifting the peak of stress due to 
temperature gradient awa^- from the rim is possible. The 
stress-concentration factor caused by the method of blade 
attachment would then apply to lower nominal stresses. 

INFLUENCE OF ELASTIC CONSTANTS 

Calculations of centrifugal and thermal stresses were made 
for several assumed variations in elastic properties to indicate 
the changes that can lie made in disk stresses by altering 
material properties and to show the effect on stress calcula- 
tions of making the simplifying assumption that the elastic 
properties do not vary with temperature. Because the 
centrifugal-stress calculations were based on a blade stress 
at the rim of 8500 pounds per square inch, this value of 
radial centrifugal stress exists at the rim for all the calcula- 
tions. In conformity with the assumed conditions, the 
radial thermal stress is zero at the rim for all variations in 
elastic properties. For a solid disk, the radial and tan- 
gential stresses are equal at the center. For all the calcula- 
tions made to investigate the effects of changes in elastic 
properties, critical conditions were found to occur only at the 
rim and at the center, and the response of the stresses to 
changes in elastic properties may therefore be observed by 
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consideration of only the tangential stresses at the center 
and at the rim. ALL calculations made to show the effects of 
change in elastic constants were based on the temperature 
distribution of figure 2 (a) with n equal to 2. 

MOnULL'S OF EL.4STiaTY 

The influence of modulus of elasticity ou centrifugal and 
thermal stresses was investigated for assumed relations 
bctwt'eu modulus of elasticity and temperature, as shown in 
figure 7. Line b approximates the actual variation of 
modulus of elasticity of a typical heat-resisting alloy in the 
t(‘niperature range shown. Centrifugal-stress calculations 
were made for the relations represented by lines a, b, and 
c of figure 7 (a) to show the effect of a change in slope of the 
modulus-temperature characteristic. The relation of the 
tangential centrifugal stress to the slope of these lines is 
shown by figure S (a). It is seen that a calculation which 
makes no allowance for the slope of the modulus of elasticity- 
temperature cmve would yield rim stresses shghtly higher 
than the correct values, but the center stresses would be 
approximately correct. 

As shown in the appendix, any calculation of centrifugal 
stresses with the assumption that the modulus of elasticity 
does not vary with temperature gives results that are inde- 
pendent of the modulus of elasticity. Calculations of 
centrifugal stresses were made for the relations represented 
by lines b, f, and g of figure 7 (c) to show the effect of a shift 
in level of the modulus-temperature characteristic at a given 
slope. The effects on stresses are shown in figure 8 (b); 



(aj Centrifugal stress for various slopes. Centrifugal stress for constant, typical 

<.c) ^.Thermal stress for various slopes. slope. 

id) Thermal stress for constant moduli. 

Figure" 8.— Effect of varying moduli of elasticity as shown in figure 7 on tangential stresses 
at center and rim of turbine disk. Speed, 11,000 rpm; blade loading at rim, 8500 pounds 
inch; center temperature, 600^ F; rim lemp^’aiure, 12X)'^ F; temperature dis- 
tribution given by r=ar’'-r&t 0 * 2 . 
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for a given slope, the center and the rim centrifugal stresses 
are negligibly affected tlie level of the modulus of 
elasticity. 

Lines a, b, and c of figure 7 (a) show assumed relations 
between elastic modulus and temperature in wliich the 
average value of the modulus is constant but the slope of 
the curve is changed. The effect on thermal stresses of such 
change is showm by figure 8 (c) ; a thermal-stress calculation, 
which assumes a constant modulus equal to the average of 
the actual moduli across the disk, would yield correct center 
stresses but the magnitude of the rim stresses would be sig- 
nificantly increased. The influence of level of the modulus 
of elasticity on thermal stresses was investigated for the 
relations represented by hues c, d, and e of figme 7 (b). The 
result is shown in figure 8 (d) where the thermal stresses are 
directly proportional to the magnitude of elastic modulus. 

COETnOENT OF THERM.YL EXPANSION 


The assumed relations between coefficient of thermal 
expansion and temperature, which were used in investigating 
the influence of this coefficient on thermal stresses, are given 
in figure 9. Line b represents the actual values of a typical 
heat-resisting alloy. Lines a, b, and c of figure 9 (a) show 
assiuned relations between coefficient and temperature in 
which the slope varies but the average value remains con- 
stant. The effect of this type of variation on the calculated 
stresses is shown in figure 10 (a), where it is seen that a 



Temperature^ ®F 

CaJ Various slopes. 

(b) Zero slope. 

(c) Constant, typical slope. 

FiGiTRE 9.— Assumed relations between coefficient of thermal eipansion and temperature 
for calculations of thermal stresses. 
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Slope, [f;>7.//n;/T]/"F 



Average of coeff/cienfs of expansion, (inlin.)/°F 


(a) Thermal stress for various siojies. 

Cb) Thermal stress for two constant slopes. 

Figuhs 10.— Effect of varying cooffldents of expansion as shown in figuj'e 9 on tangential 
stress at center and at rim of turbine disk. Speed, 11,000 rpm; blade loading at rim, 8500 
pounds per .square inch; center temperature. 600° F; rim temperature. 1200° F; temperature 
distribution given by T=ar»-H>.' 


calculation of thermal stress, which neglects the variation of 
coefficient of expansion with temperature, yields stresses 
slightly lownr than the correct values. The sets of the values 
represented by lines c, d, and e of figure 9 (b) were used to 
determine the effect of magnitude of the coefficient when it 
does not vary with temperature. Figure 10 (b) shows that 
the thermal stresses are directly proportional to the level of 
coefficient. The sets of values represented b}" lines b, f, and 
g of figure 9 (c) were used to determine the effect of magni- 
tude of coefficient of expansion when the coefficient of 
expansion-temperature characteristics have a given slope. 
Figure 10 (b) shows that a linear relation exists between the 
level of the coefficient and the stresses. 

POISSON’S RATIO 

Calculations of thermal and centrifugal stresses were also 
made for values of Poisson's ratio of 0.2, 0.3, and 0.5. Most 
heat-resisting alloys have a value of Poisson’s ratio at room 
temperature of approximately 0.3. At elevated tempera- 
tures, a material approaches the plastic condition and, ac- 


cording to the theory of plasticity, Poisson’s ratio for the 
fuUy plastic state would be 0.5 (reference 3, p. 79). Thus, 
the range of 0.2 to 0.5 covers all values of Poisson’s ratio 
that may be ex'pected to occur. 

A change in the value of Poisson’s ratio from 0.2 to 0.5 
decreased the centrifugal tangential stress at the rim by 22 
percent and the thermal stress at the rim by 1.4 percent and 
increased the net or total stress at the rim by 6 percent. 
For disks in which thermal stresses are predominant, the 
precise evaluation of Poisson’s ratio is therefore unimporlant. 

SUMMARY OF RESULTS 

When perfect elasticity was assumed, the following quali- 
tative results concerning the influence of temperature gra- 
dients and elastic properties on stress conditions in aircraft 
gas-tufhuiG disks were obtained: 

1. Severe temperature gradients caused thermal stresses 
of sufficient magnitude to reduce the operating safety of the 
disk. Attempts to increase the strength of the material b}' 
cooling the center of a disk might reduce the margin of 
safety as a result of the thermal stresses set up. 

2. Small temperature gradients were desirable because 
they created thermal stresses that subtracted from the 
centrifugal stresses in the region of the rim. Center stresses 
were increased by the small temperature gradients. 

3. Thermal gradients produced a tendency for the maxi- 
mum stress to exist near the rim. By altering the tempera- 
ture distribution, the position of this maximum stress could 
be shifted away from the rim. A stress-concentration factor 
caused by the method of blade attachment then applied to 
lower nominal stresses. 

4. Centrifugal stresses were slightly aiTected by changes 
in the relation between modulus of elasticity and tempera- 
ture. The thermal stresses were approximately proportional 
to the coefficient of expansion and to the modulus of elas- 
ticity. Because the olTocts of thermal stresses arc likt'ly to 
predominate, iPis desirable that the modulus of elasticity 
and coefficient of thermal expansion be small. Design cal- 
culations require accurate data on clastic modulus and coef- 
ficient of thermal expansion over the temperature ranges 
encountered. 

5. A change in Poisson’s ratio from 0.2 to 0.5 decreased 
the centrifugal tangential stresses at the rim but increased 
the total stresses bj^ a negligible amount. Elevated- 
temperature measurements of Poisson’s ratio were therefore 
considered to be of small value for calculations of gas-turbine- 
disk stresses. 


Flight Propulsion Research Labor.^tory, 
National Advisory Committee for Aeronautics, 
Cleveland, Ohio, March 29, 1947. 





APPENDIX 


MATHEMATICAL ANALYSES 


SYMBOLS 

The following S 3 anbols are used in this analysis: 

-1 maximum shear-strain enei-gy, in.-lb/cu in. 

E modulus of elasticity, Ib/sq in. 

G modulus of rigidity, Ib/sq in. 

h axial thickness of disk, in. 

r disk radius, in. 

,v equivalent tensde stress, Ib/sq in. 

«ni stress margin of safety, Ib/sq in. 

.S'n tensile test stress, Ib/sq in. 

Si, S 3 principal stresses, Ib/sq in. 

AT difference between temperature of disk and tem- 
penature at zero thermal stress, °F 
Cl mean coefScient of thermal expansion, (in./in.)/°F 

11 Poisson’s ratio, 0.3 

p mass density of disk material, lb secVm.‘‘ 

0 -, total radial stress, Ib/sq in. 

at total tangential stress, Ib/sq in. 

u angular velocity, radians/sec 

EQUIVALE.NT TEXSnE STRESS 

The maximum shear-strain energy is considered to be the 
index of the severity of the stress state in a high-temperature 
rotating disk. This quantity is given in reference 3 (p. 73) 
as 


= [(Si — 82)H(S2-S3)^H-(S3 -Si)^] (1) 


In a rotating disk, (Xr and are principal stresses whereas the 
axial stress is zero, so that at any point 

A.=^^[{a-r—a-,y^-\-at^+<Tr‘] ( 2 ) 

The stress s in a tensile specimen is rmiaxial; the substitution 
in equation (1) of s=Si and S 2 =Sa =0 gives 


( 3 ) 


The uniaxial tensile stress, which is equivalent to the biaxial 
stress in a rotating disk, is now defined by equating these 
two values of maximum shear-strain energy. Thus 

[(cr— 0-j^+ a-r°] (4) 


or 


S= x/a-r^— crrO-t+ crt C5) 

Then, inasmuch as is the stress-rupture value of a disk 

material at a temperature corresponding to the temperature 
in a disk at a point where the state of stress is being eval- 
uated, the margin of safety 8„ is defined as 

Sm=8o — s (6) 

EFFECT OF CONSTANT MODULUS OF ELASTiaTY ON CENTRIFUGAL 

STRESSES 

The equilibrium equation (7) and the compatibility equa- 
tion (8) are given in reference 1. 

|;(rAa,)-Aff,+pco*Ar2=0 (7) 


d (St_ d ixctt , d 


drETr E^dr ^ 


(1-h/x) (g-f— 0-,) 

Er 


( 8 ) 


Equation (7) contains no terms that involve the modulus of 
elasticity. In equation (8), the substitution of a=0 and a 
constant value of modulus of elasticity establishes that the 
centrifugal stresses are not influenced by the magnitude of 
any constant modulus of elasticity because it cancels out of 
the compatibility equation and is not involved either in the 
equilibrium equation or in the boundary conditions. 
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